Reponse to Comment on "Direct mapping of the finite temperature phase diagram of 

strongly correlated quantum models" 


Qi Zhou 1 , Yasuyuki Kato 2 , Naoki Kawashima 2 , and Nandini Trivedi 3 

1 Joint Quantum Institute, University of Maryland, College Park, MD 20742 
2 Institute for Solid State Physics, University of Tokyo, Kashiwa, Chiba 277-8581, Japan 
3 Department of Physics, The Ohio State University, Columbus, OH 43210 


Pollet et al have questioned our claim staked in pQ that 
from the location of sharp features in the compressibility 
k in the trap, it is possible to map out the phase bound- 
ary between two phases within reasonable accuracy. We 
stand by it and use this reply to sharpen the conditions 
under which "sharp" features in k reflect the effects of 
critical fluctuations. 

(1) The first point raised by Pollet et al [2 is that 
sharp density gradients may also lead to sharp features 
in n. We agree both critical fluctuations and strong den- 
sity gradients, in general, can contribute to the features 
in ft, as seen in Fig. 4 of our PRL. However, it is easy 
to separate out effects due to density gradients by simply 
making the trap potential flatter or by making the system 
larger by increasing the number of particles. The effects 
due to density gradients will then be diminished. It can 
also be checked theoretically, by closely corroborating the 
kink features, with the superfluid density distribution us- 
ing the local density approximation, as was done in Fig. 5 
of our PRL. 

(2) The second point raised by Pollet et al is that when 
they simulate the parameters in our Fig.(5B) in a trap, 
they observed no features in k associated with the critical 
point by showing n within a small region of about 7 lat- 
tice spacings near the critical point r c in their comment. 
We do not dispute their numerical data but we disagree 
with their analysis and interpretation. 
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We are well aware that singular features in k will get 
rounded by finite size effects. However, we want to make 
a definite point here, that even with this rounding, if the 
data is analyzed carefully, it should be possible to locate 
the phase boundary within reasonable accuracy. 

We first demonstrate this fact using the same data gen- 
erated by Pollet et al. In Fig. 1(a), We re-plot their data 
in a larger region, up to 17 lattice spacings, than what 
has been shown in their comment. It is now clear that 
there is a maximum in k from which r c can be easily lo- 
cated. On a larger scale shown in Fig. 1(b), the structure 
of n becomes even clearer. 

Next we address the question about the accuracy of 
the method for extracting the phase boundary ji c from 
the maximum in n. A quantitative determination of \i c 
and the associated error bar can only be obtained from a 
detailed finite size scaling analysis, however, a qualitative 
estimate can be made by the procedure outlined below. 
In Fig. l(c-d), we show our quantum Monte Carlo simula- 
tions in the trap. The finite size rounding in a trap gives 
an uncertainty of r c /d between 29.6 and 30.4, where d is 
the lattice spacing. This leads to an uncertainty in the 
critical chemical potential /jl c /U = —0.094 ± 0.016. Our 
result for the mean value of ji c is within 2% compared 
with the mean value of ji c obtained in [3] at the same T/t 
and t/U. The critical density n c located as 0.375 ±0.035 
is also within 1.4% of the mean value of n c found in [3]. 

Thus based on this comparison with Pollet et al's own 
results, we emphasize that the finite temperature criti- 
cal point of quantum models can be identified within a 
reasonable accuracy from the derivative of the density 
profile of atoms in the trap, despite the finite size round- 
ing effect in the vicinity of the critical point. 

We also use this opportunity to correct typos in our 
PRL pointed out by Pollet et. al: The x label of Fig.(4A) 
should be 10, 20, 30, 40 and also the triangle and box sym- 
bols should be reversed. The label of fi/U for Fig.(5B) 
should be 0.5, 0.448, 0.290, 0.028, -0.339. 


FIG. 1: (a-b) Re-plot, from Fig. 1 of the comment [2], of the 
trap density profile n and local compressibility k over a larger 
scale. (c-d)Our direct simulations in a trap for the same T 
and t/U as in (a-b) but at different trapping frequency V/t — 
0.00262 and chemical potential fi/t = 1.984. 
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